home *** CD-ROM | disk | FTP | other *** search
/ Enter 2004 January / enter-2004-01.iso / files / maxima-5.9.0.exe / {app} / share / maxima / 5.9.0 / demo / eaton2.dem < prev    next >
Encoding:
Text File  |  2003-02-09  |  804 b   |  27 lines
  1. SHOWTIME:TRUE$
  2. LOADPRINT:FALSE$
  3. EQN1:A*COS(P2)=S+B*COS(P3);
  4. EQN2:A*SIN(P2)=E+B*SIN(P3);
  5. /* Declare p2, p3 and s to be time-dependent. */
  6. DEPENDS([P2,P3,S],T);
  7. /* [1] From equations 1 and 2, eliminate P3: */
  8. EQN3:EXPAND((EQN1-S)^2+(EQN2-E)^2);
  9. EQN3:TRIGSIMP(EQN3);
  10. /* [2] Using the results of step 1, solve for S in terms of
  11.        a,b,e,p2: */
  12. S_SOLUTION:SOLVE(EQN3,S);
  13. /* Note that this differs from the result given.
  14.    [3] Take the derivative of eqn2 wrt t:  */
  15. EQN3:DIFF(EQN2,T);
  16. /* [4] Solve for p3-dot in terms of a, b, p2, p2-dot,p3. */
  17. P3_DOT:SOLVE(%,DIFF(P3,T));
  18. /* [5] Take the derivative of eqn2 wrt t:  */
  19. EQN5:DIFF(EQN1,T);
  20. /* [6] Solve for s-dot in terms of a, b, p2, p2-dot,p3. */
  21. SOLVE(EQN5,DIFF(s,T));
  22. %,P3_DOT,FACTOR;
  23. TRIGREDUCE(%);
  24. /* [7] */
  25. Z:I*(A*COS(P2)/(B*COS(P3)))^2;
  26. DIFF(Z,P2);
  27.